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Abstract. In this paper we study a model describing the growth of 
necrotic tumors in different regimes of vascularisation. The tumor con- 
sists of a necrotic core of death cells and a surrounding nonnecrotic 
shell. The corresponding mathematical formulation is a moving bound- 
ary problem where both boundaries delimiting the nonnecrotic shell are 
allowed to evolve in time. We determine all radially symmetric station- 
ary solutions of the problem and reduce the moving boundary problem 
P^ ' into a nonlinear evolution. Parabolic theory provides us the perfect con- 

•^D . text in order to show local well-posed of the problem for small initial 

* ' data. 
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1. The mathematical model 
>. 
*sO In this paper we study a moving bondary problem describing the growth 

j*^ ' of a necrotic tumor in the absence of inhibitor. The model purposed ini- 

^Nl ■ tially in [4,15,17], was reformulated by using algebraic manipulations [5,12] 

!— " | to describe evolution of tumors in all regimes of vascularisation. Neverthe- 

^^ ■ less, the analysis in [5, 12] is simplified by the assumption that the tumor 

core is nonnecrotic. This aspect is considered into our modeling, where fol- 
lowing [10, 17, 18], we assume that the tumor consists of a core of death 
cells (necrotic core) and a shell of life-proliferating cells surrounding the 
core (nonnecrotic shell). The blood supply provides the nonnecrotic region 
!-h ' with nutrients, while there is no blood supply in the necrotic region and 

the concentration of nutrients is at a constant level which cannot sustain 
cell proliferation. However the model presented here includes two moving 
boundaries, one parametrising the boundary of the necrotic core and one 
for the outer boundary of the tumor, both of them having infinitely many 
degrees of freedom. This fact makes the problem more involved in com- 
parison to other models which either neglect the necrotic core [3, 6-10] or 
consider only the radially symmetric problem when the tumors are annular 
domains [10,17,18]. 



>< 



Ctf 



2000 Mathematics Subject Classification. 35K55; 35R35; 35R37. 

Key words and phrases. Radially symmetric; Stationary solution; Classical solution; 
Fourier multiplier. 

1 



J. ESCHER, A.-V. MATIOC, AND B.-V. MATIOC 
The mathematical model is given by the following system of equations 



' Aijj 


= 


t 




in tt(t), 


£> 


Ap 


= 







in Q(t), 


t> 


Y> 


= 


G 




on Ti(t), 


t > 


Y> 


= 


G-Vo 




on r 2 (t), 


t > 


p 


= 


\x\ 2 
«ri(t) -^G— 




on Ti(t), 


t > 


p 


= 


|x| 2 
K r 2 (t) - AG — ~ 


V'o 


on r 2 (i), 


i > 


Vi(t) 


= 


d v ^-8 Vi p-AG 


Vi ■ X 
2 


on Ti(t), 


t > 


, n(o) 


= 


Qq, 








re fl(t) 


Z 1 


2 is the domain occupied by the nonnecro 



:i.d 



1,2 



rate at which nutrient is added to Q(i), over the outer boundary Ti(t), by 
the vascularisation, p is the pressure, r 2 (t) is the interior boundary enclosing 
the necrotic core, V{ the outward orientated normal and «r,(t) the curvature 
of Ti(t), i = 1, 2. By convention, ftn(t) is positive and n^ 2 t t \ negative if r^i), 
1 < i < 2 are close to a circle. Moreover, VJ(t) stands for the normal velocity 
of Ti(t), the constants A, G € M have biological relevancy being related to 
cell proliferation, cell apoptosis, and vascularisation. The scalar ipQ > is 
linked with the constant nutrient concentration assumed within the necrotic 
region. The initial tumor domain is given by Qq. 

For a precise deduction of the system (1.1) and its biological meaning we 
refer to [5, 12], the only difference to the model presented there being the 
consideration of the interior necrotic region bounded by r 2 (i). 

The first main result of this paper is the following theorem: 

Theorem 1.1. Given (i?i,i? 2 ) € (0, oo) 2 with i? 2 < R\, let ipQ be the constant 
defined by (2.8). There exists A £M and GgR\ {0}, such that the annulus 

A{R l ,R 2 ) := {x G R 2 : R 2 < \x\ < Rt}, 

is a stationary solution of problem (1.1) provided ipo ^ i/jq. Moreover, A and 
G are uniquelly determined by i?i,i? 2 , and i/jq. 

If G = 0, then problem (1.1) has no radially symmetric stationary solu- 
tions. 

In contrast to [5,12], where the radially symmetric stationary tumors are 
circles with radius which depends only on the constant A, the radii of the 
stationary annular tumors found in Theorem 1.1 depend on both constants 
A and G, but also on i/jq, cf. (2.9). 

In order to prove local-well-posedness of the moving boundary problem 
(1.1) the well-posedness of system (1.1) (see Theorem 1.2 below) we introduce 
first a parametrisation for the interfaces Ti(t) and T 2 (t), which are the main 
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FIGURE 1 . Parametrisation of the tumor domain 



unknowns of system (1.1). Let < R 2 < R\ be given and fix a G (0, 1). We 
set 

V := {p G h 4+a (S) : \\p\\ C ($)<a}, 

where 

R\ — R 2 

a < . 

R\ + R 2 

The small Holder space /i m+/3 (§), f3 G (0,1) and m G N, is defined as the 
completion of the smooth functions in C m+ ^(S). Each pair (pi, P2) G V 2 
parametrises a C 4+Q -domain 

n( Pl ,p 2 ) := {y G R 2 : R 2 (l + p 2 (y/\y\)) < \y\ < i?i(l + pi (y/\y\))} . 

The condition on a ensures that the boundary portions of 0(px,P2) 

r(pi) := {x : \x\ = Ri(l + pi(x/\x\))}, 

i = 1,2, are disjoint (see Figure 1) for any choice of (pi,p 2 ) G V 2 . They can 
be seen to be zero level sets, T(pi) = i\T pi — 1 (0), where i\T w : M 2 \ {0} ->• R, i = 
1,2, are defined by 

N Pi (x) = \x\ - Ri - RiPi (x/\x\) , x / 0. 

Hence, the outward unit normal at dQ(pi,p 2 ) is given by 
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If (pi,p 2 ) '■ [0,T] —7- V 2 describe the motion of the tumor, then we can 
express the normal velocity of both boundary components in terms of pi by 
the formula 
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With this notation, system (1.1) becomes a problem with p-y and p 2 as 
unknowns: 
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(1.2) 
A pair (p\, p2,i^,p) is called classical solution of (1.1) on [0,T],T > 0, if 

Pi € C([0,T],V)nC\[0,T],h 1+a (S)), i = 1,2, 

ifj(t,-),p(t,-)£buc 2+a (n(p 1 (t),p 2 (t))),te[0,T], 

and if (pi,p2,ifi,p) solves (1.2) pointwise. Given U C M 2 open, we set 
6wc +a (C7) to be the closure of the the smooth functions with bounded and 
uniformly continuous derivatives BUC°°{U) within BUC 2+a (U) (if U is also 
bounded then BUC 2+a {U) = C 2+a (U).) 

Concerning well-posedness of system (1.1), our second main result states 
that problem (1.1) possesses a unique solution provided that initially the 
tumor is close to an annulus (which is not necessarily a stationary solution). 

Theorem 1.2 (Local well-posedness). Let < R 2 < R\, and {A,G,iJ)q) G M 3 
be given. 

There exists an open neighbourhood O C V such that for all (pi,p 2 ) € O 2 , 
problem (1.2) possesses a unique classical solution defined on a maximal 
time interval [0, T(poi) P02)) an ^ which satisfies (pi, p 2 )(t) S O 2 for all t G 
[0,T(poi,Po2)). 

The outline of this paper is as follows: we study in Section 2 the radially 
symmetric free boundary problem which describes the stationary solutions 
of problem (1.1) and prove Theorem 1.1. In the last section of this paper we 
prove the local well-posedness result Theorem 1.2. 

2. Radially symmetric stationary solutions 

We determine in this section the radially symmetric steady-state solutions 
of (1.1), situation when the nonnecrotic shell is a steady annulus. 

The most simple situation is the case G = 0, when the problem is invariant 
under translation and rotations.. Then, the annulus A(R\,R 2 ) centred in 
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zero with radii Ri > R 2 , is a stationary solution of system (1.1) if and only 
if 

P '(R l ) = iP'(R i ), i = l,2, 
where p is the solution of the problem 

p" + -p' = 0, R 2 <r<R u 

pik) = Ri~ l - AGRl/4, (2-1) 

p{R 2 ) = -R 2 ~ 1 -AGR 2 2 /A-^ 

when G = 0. System (2.1) corresponds to the Dirichlet problem for the 
pressure p in (1.1) (the second, fifth, and sixth equations of (1.1)), where 
we used polar coordinates when expressing the Laplacian. Notice that the 
boundary data are constants, thus p depends only on r, the distance to the 
origin. 

Given GeR, the solution of (2.1) is given by the relation p(r) = a RlR2 ln(r)+ 
&R1R2) r 2 < l r l < Ri, w ith 

^r 1 + R2' 1 + AG (Rj - Rf) /4 + jjp 
aRlR2 ~ HR1/R2) 

b RlR2 = Rr 1 - AGR 2 JA - a RlR2 ln(R 1 ). 

Furthermore, ip is the solution of the problem 

ip" + -1// - ij; = 0, R 2 <r<R x , 

7 ' Tp(Rx) = G, (2-2) 

ij(R 2 ) = G-Vo, 

when G = 0. Also, for fixed G E R, the solution of (2.2) can be written 
as linear combination of modified Bessel functions of first and second kind 
V> = ckft/o + c 2 RlR2 K , with scalars 

i GK (R 2 ) + (tf; -G)K (R 1 ) 



' RlRa I (Ri)K (R 2 ) - I (R 2 )K (R 1 ) ' 

2 -GJ CR2)-M)-G9lo(fli) 

CRlR2 I (Ri)K (R 2 ) - 10(^2)^0(^1) ' 

Consequently, A(R\, R 2 ) is a steady-state solution of (1.1) when G = if 
and only if 

— — 

fii + fi 2 + ^ l __,, KojR^hjRi) + I (R 1 )K 1 (R l ) i = 12 (23) 

ln(Ri/R 2 ) Ri ^ I (Ri)K (R 2 ) - I (R 2 )K (RiY >>^-> 

where we used the relations / = I\ and K' = —K\. It follows then easily 
that the system consisting of the equations (2.3) has solutions (R\,R 2 ) with 
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R\ > R 2 exactly when 

R2 ir (fli)Ji(fli) + Jo(fli)i<ri(fli) (2 4 

R ± K (Ri)h(R2) + I {Ri)K 1 {R 2 y [ ~> 

Equation (2.4) is obtained by expressing ipQ in both relations (2.3) and setting 
them to be equal. We show now that equality holds in the relation above 
only when R\ = R 2 . Indeed, fix R± > and consider the auxiliary function 
g : (0,Ri\ ^R with 

g{x) = K Q {R 1 )xh{x) + h{R 1 )xK 1 {x) - R 1 (K (R 1 )I 1 (R 1 ) + I (Ri)K 1 (R 1 )) 

for < x < R±. Obviously g{R\) = 0. If we show that the derivative g' has 
constant sign on (0, R\] then we are done, that is there is no positive R 2 < R\ 
such that (Ri, R 2 ) solves (2.4). Well-known properties of the modified Bessel 
functions (see [2]) lead to 

g'(x) =K (Ri)h(x) + I (R 1 )K 1 (x) 

+ K (R 1 )x(I (x) - {l/x)h{x))+I Q {Rx)x{-K Q {x) - (l/x)K 1 (x)) 
=x{I Q (x)K {R l ) - I Q (R{)K Q {x)) < 

for all x € (0,R\). That the last expression is negative is a consequence of 
the following facts: Iq and Kq are both positive functions, Iq is strictly in- 
creasing, and Kq is strictly decreasing. Hence, problem (1.1) has no radially 
symmetric stationary when G = 0. 

Let now G 7^ 0. In this case A(R\,R 2 ) is a steady-state solution of (1.1) 
exactly when 

ip RlR2 (Ri)-p'(Ri)-AG^ = 0, i = l,2. (2.5) 

Using again the relations I Q = I\ and K Q = —K±, the identities (2.5) re- write 

c RlR M R i) ~ 4l lR2 Ki(Ri) - URiR*^. ~ AG^ = 0, i = 1,2, 

which seem to be very involved as expressions of variables Ri and R 2 when 
trying to solve the system consisting of both of them. However, they can be 
viewed as equations for A and G 

aiG + biAG = Ci, i = 1,2, (2.6) 

with coefficients aj,6j, and q given by: 

(K (R 2 ) - KoiR^hiRi) - (Jo(fli) - I (R 2 ))Ki(R l ) 



a, ; : = 



I (Ri)K (R 2 ) - I (R 2 )K (R X ] 

R^ — i?2 1 Ri 



bi ''~ Ahx{R 1 /R 2 )R i " 2 ' 



. KoiR^hiRi) + Ioi^KiiRi) , R 1 - 1 + R 2 - 1 +^ 1 

Ci := -tpo = , „ > „ -, = ; r - n r- T —- rFr T + 



' I (Ri)K (R 2 ) - IoiR^KoiRi) In(i2i/i2 2 ) R t 
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Figure 2. The expression a\b 2 — «2&i and ipQ, for fixed R±, 
as a function of the variable R 2 £ (0, R\). 



The system of equations (2.6) has a (unique) solution (A, G) with G /0 
provided that 



ai& 2 - a 2 bi / 0, ci6 2 - c 2 b\ ^ 0, 

and c\ 7^ or C2 7^ 0. 



(2-7) 



The computation done for the case G = shows that ci and c 2 cannot be 
simultaneously zero when R 2 < R\ . For fixed R\ > we may see the expres- 
sion a\b 2 — a 2 b\ as a function of R 2 £ (0,i?i). This function is decreasing 
with respect to R 2 (see Figure 2), thus 0162 — o, 2 b\ only when i?i = R 2 . 
Furthermore, b\c 2 = b 2 c\ if and only if ^0 = ^o; where 



V'n 



( bl / Rl - b 2 /R 2 ) ll*0§ 



K (R 1 )(b 1 I 1 (R2)-b 2 I 1 (R 1 ))+I (R 1 )(b 1 K 1 (R 2 )-b 2 K 1 (R 1 )) 
Io(R 1 )K (R 2 )-I (R 2 )Ko(Ri) 



+ 



Rf-3 



2R 1 R 2 ln(R 1 /R 2 ) 

(2.8) 
It is not difficult to see that the numerator of the fraction is negative, and 
the same holds true for the denominator, implying that tpQ > 0. We plotted 
in Figure 2 the expression on the right hand side of (2.8) for fixed R\ > in 
dependence of R 2 € (0, R\). Consequently, A(R\, R 2 ) is a stationary solution 
of (1.1) if and only if tpQ is not the critical constant given by (2.8) and 



a\c 2 — a 2 c\ 



c\b 2 
This proves Theorem 1.1. 



C2&1 



G 



c\b 2 - c 2 bi 
a\b 2 - a 2 b\ 



(2.9) 



3. The moving boundary problem 

This last section is dedicated entirely to the proof of our secon main re- 
sult, Theorem 1.2. In order to prove well-posedness of problem (1.1) in the 
context defined in the introduction we transform first (1.2) into a problem 
on the fixed domain Q := (7(0,0), with boundary Ti := i?i§ and T2 := R 2 S>. 
This transformation will allow us to introduce solution operators related 
to problem (1.2) and which will enable us to reduce system (1.2) into an 
abstract nonlinear evolution equation for the pair (pi,p 2 ). 
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Pick therefore < R 2 < Ri, (A,G,ip ) G M 3 , and a G (0,1). Given 
{pi, P2) € V 2 , we define the mapping ® Pl , P2 '■ fi — > ^(pi,p 2 ) by the relation 

{Rx - \x\)R 2 (l + p 2 (x/|x|)) + (|x| - R 2 )Ri(l + Pi(x/\x\)) x 



"pijisW 



i?i — i?o 



for 2; G £1. One can easily check that Q pitP2 is an diffeomorphism, i.e. 
®Pi,P2 ^ Diff 4+a (Q, Q(pi,p 2 )), which maps Tj onto T(pi),i = 1,2. Using 
this diffeomorphism, we define the transformed operators 

A(pi,p 2 ) : buc 2+a (n) -> 6 uc a (0), >l(pi, p 2 )« := A (v o 9^ p2 ) o9 pi , P2 -v, 

which is an elliptic operator depending analytically on (pi,p 2 ), i.e. 

A g c u (v 2 , c{buc 2+a (n), buc a (n))), (3.1) 

respectively the trace operators B{ : V 2 x (buc 2+a (Q)) 2 — >■ /i 1+a (§) by 

Bi(pi,p2,v,q) := —Ci(pi,p 2 )v- — Ci(pi,p 2 )q- Vi{p u p 2 ). 
Hi Hi 

Given (pi ,/) 2 )eV 2 , the linear operators C i (p 1 ,p 2 )e£( buc 2+a (0) ) , /i 1+a (S) ) , 
2 = 1,2, are given by 

Ciipupdviv) := (V^oe-^l ViV Pi )oe pi)P2 (i?,y) 

for v G 6uc +a (17) and y G S. Moreover, 

A(pi,p 2 ) -=-^-(2 VAr ft)°@Pi,p 2 (^y)- 

The operators Cj and Pj,i = 1,2, depend analytically on (pi,p 2 ) too, 

d G C W (V 2 , £(?mc 2+a (ft)), /i 1+a (§))) and V { G C7 W (V 2 , /i 1+Q (S)). (3.2) 

Having defined this operators we may re-write now (1.2) in an equivalent 
form. Namely, if (pi, p 2 ,ip,p) is a solution of (1.2), v := ip o @ pltP2 , and 
q := p o pi:P2 , then the tupel (pi, p 2 ,v, q) solves the following system 



A(pi,p 2 )v 


= 







in 0, 


t>0, 




A(pi,p 2 )q 
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in 0, 


t>0, 
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on Ti, 


t>0, 
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G-i> 




on r 2 , 


£>0, 
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^(Pi)-^(1 + Pi) 2 




on Ti, 


t >0, 
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-i^ 2 )-^(l + p 2 ) 2 - 


-tpo 


on r 2 , 


t>0, 




dtPi 
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&i(pi,P2,v,q) 




on §, 


t>0,i = 


= 1,2 


Pi(0) 


= 


P01, 










P2(0) 


= 


P02, 











(3.3) 
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where k : V — > /i 2+ "(§) is defined by 

, * (l + p) 2 + 2p a -(l + pK 

* (P) - = ((l+p)2 +/ y2)3/2 ' " GV > 

and we identified functions on T{ with those on §, i = 1,2, via the difieo- 
morphisms [§ 9 y h-> i^y G Tj]. 

Though the problem becomes more involved (the diffeomorphism intro- 
duces additionally nonlinearities), (3.3) has the advantage that the sets where 
the differential equations and the boundary conditions are imposed do not 
change with time. It is convenient now to introduce solution operators to 
Dirichlet problem closely related to system (3.3). 

Lemma 3.1. Given (p\,p2) £ V 2 , we let T(pi,p2), S(pi,p2) € buc 2+a (Q) 
denote the unique solution of 

A(pi,p2)v = v in Q, 

v = G onTi, (3.4) 

v = G-V>o on T 2 , 

and 

A(p\,p2)q = in O, 

q = ^{pi)-^{l + Pif onV u (3.5) 

Q = ~^«(P2) 4^(! + P2) 2 -i>o onT 2 , 

respectively. The operators T and S depend analytically on (pi,p2)- 

Proof. Given {p\,p2) € V 2 , the Dirichlet problems (3.4) and (3.5) are uniquelly 
solvable, cf. [16, Theorem 6.14]. Moreover, since A and k depend analyti- 
cally on their variables we deduce that also T and S do that. We may take 
now into consideration that T and S map both the smooth functions into 
BUCr +a (£l) and conclude that their range is contained in buc +a (U). □ 

With this definion (3.3) reduces to the following evolution equation 

8tX = $(X) X(0) = X , (3.6) 

where X := (pi,p2), Xq := (poi,po2), and $ := ($i,<J>2)- The components 
of the nonlocal and nonlinear operator $ are defined as follows 

$i(pi,P2) ■■= Bi(Pi,T(pi,p 2 ),S(pi,P2)), i = 1,2. 
In order to prove well-posedness of problem (3.6) it suffices to show that 

a Pl $i(o) d P2 ^(o) 

d Pl <P 2 (0) d P2 <P 2 (0) 



89(0) 



generates a strongly continuous and analytic semigroup. The key role is 
played by the operator S which depends on the highest order derivatives of 
Pi,i = 1, 2. We have: 
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Theorem 3.2. The operator $ is analytic, i.e. <£ G C W (V 2 , (h 1+a (§)) 2 ) . Green 
/? £ (0, 1), the Frechet derivative <9<I>(0), seen as an unbounded operator in 
(/i 1+ ^(§)) 2 with domain (/i 4+ "(§>)) 2 generates a strongly continuous and an- 
alytic semigroup in £((/i 1+ ^(§)) 2 ), i.e. 

-d$(0) G ^((/i 4+/3 (S)) 2 , (/i 1+/3 (S)) 2 ). 

Proof. The regularity assumption follows directly from (3.1) and (3.2). More- 
over, since the constant a fixed at the begining of this section was arbitrary, 
we may replace a by (3 and all the assetions already established remain valid. 

Let us now study the Frechet derivative of $ in 0. One can easily see that 
the highest order terms in (pi,p2) of d&(0,0)[(pi, P2)] are those obtained 
when differentiating the curvature operator. 

Consider first <9 pi $i(0). Since 3/«(0)[p] = -p" - p for p G p 4+a (§), we 
may decompose 

d Pl *i(0)[f>i]=A 11 + B 11 , 

where B\\ is an operator of first order, i.e. B\\ G C{h 2+ "(B>), /i 1+ ^(§>)), 
1 



Anpi := -- f C 1 (0)(A,tr 1 ,tr 2 )- 1 (0, /) ' 1 ',0), V Pl G /i 4+/3 (§), 



and trj,i = 1,2, is the trace operator with respect to Tj. We determine 
now a Fourier expansion for the highest order term of d pi $i(0)[hi\. Given 
p\ G /9 4+ "(§), the function w := (A, tri, tr2) _1 (0, h'{,0) is the solution of the 
Dirichlet problem 



(3.7) 



My) = ^2pi( m )y m and w(n/) = 5Z w ™( r )y k 

m m£Z 

for y G § and R2 < r < R\, we find that u>o = 0, and u; m solves, for \m\ > 1, 
the problem 





' Aw = 


in f2, 


< 


u; = p'{ 


on Ti, 




w = 


on T2. 


If we expand 







..// , 1/ m 



J m ' w m 2 



R 2 <r < R t 

Wm(Ri) = -m 2 pi(m) 
w m {R2) = 0. 



Hence 



JX2 Jt^ — JI2 -'H 



w m (r) = — — m pi(m 
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and therewith 
Aupi{y) = --^d(0)u;(y) = —^(Vw(R iy )\y) = --1— (w(ry))\ r=Ri 

1 r>\ m \r>-\ m \ , n~\ m \ pH 

= -t^ > I,, i i 7 i \m \ pi(m)y m . 

1 mGZ\{0} n l U 2 n l n 2 

We proceed similarly and write c?p 2 $i(0) = A12 + Si2> where B12 G 
£(h 2+ P(S),h 1+ P(S)) and 

A12P2 := -^=-Ci(0)(A,tr 1 ,tr 3 )- 1 (0,0,^ / ) Vp 2 G fr 4+/3 (S). 
K1JX2 

Given pi G /i 4+ "(§), the function it! := (A, tri, tr2) _1 (0, 0, p 2 ') i s the solution 
of linear Dirichlet problem 

( Aw = in Q, 

w = on Ti, (3.8) 

w = p 2 ' on T2- 

A Fourier series ansatz as we did before yields that 

T>m„—m jy—mrn 
1 \ V^ •'•'I ' - n '1 ' 9^ 1 \ m 

wiry) = — > m po(m)y 

m€Z\{0} 1 2 12 

for all y G § and R2 < r < R\, provided that /92 = /CmeZ /^(t^)?/" 1 - Whence, 
Au £ frM» m = -/^fe S ff H ff -H 2 ff -H g H H 3 feMl/ m - 

meZ ! 2 m.eZ\{0} K l ^2 ~ -^1 ^2 

We consider now the second component $2 and continue our computa- 
tion following the same scheme. The second diagonal element of the ma- 
trix 5<fr(0) may be also written as the sum d p2 &2(0) = ^22 + -E>22 5 with 
B 2 2e£(h 2+l3 (S),h 1+ P(§)) and 

A22P2 := -ic 2 (0)(A,ta lj tr 2 )- 1 (0,0, / i>2) Vp 2 G /i 4+/3 (S). 
it 2 

Using once more the expansion for the solution of (3.8), we find out that 

1 r>\ m \r>-\ m \ , p-HpH 

A22 T Pi{m)y m = — ^ V r 1 11 11 i |m| 3 ^(m)j/ m . 

zz / 4 r^\ /a „3 / j | m | _ | m | r>—\m\ pM 

m6Z 2 mGZ\{0} ""T ^2 ~ ^1 ^2 

for all p 2 = E me zP2(m) ? / m within /^(S). Finally, <9 Pl $ 2 (0) = A 21 + B 21 , 
where B 21 G £(/i 2+/3 (§), /i 1+/3 (S)) and 

4upi := -^-C2(0)(A,tT 1 ,tr 2 y 1 (0, p'1,0) Vp 2 G /i 4+/3 (S). 
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Since (A,tri, tr2) _1 (0, p'{,0) is the solution of (3.7), we may use the expan- 
sion found at that point of the proof and obtain that 

^E^Hf = -i E ff H ff -H 2 ff -H ff H H 3 Pi(m)y m 

m&L 2 meZ\{0} a l n 2 -"l n 2 

for all functions pi = ^mez/M m )y m i Q /i 4+/3 (§)- 

Let us notice that the operators A^j,l < i,j < 2, found above are all 
Fourier multipliers, since they are of the form 



J>i(m)y m .-► £ M fc j5i(m)y m 

mgZ mgZ 

with symbol (M^kez C C Using [14, Theorem 3.4], which is a theorem char- 
acterising multiplier operators between Holder spaces by studying some gen- 
eralised Marcinkiewicz conditions for the symbol of the operator, we find out 
that -An G ^(/i 4+/3 (S),/i 1+/3 (S)),i = 1,2, and that A 12 ,A 21 G £(/i 2+/3 (§)). 
This may be seen form the following relations 

9 fi>H r-M i p~M pH 

f . n aTnr | Jtj Jt 2 + -"4 ^2 . -. 

D H D -|m| „-|m| D |m| ~~Hm|->°o u auu , | ,, m | |m | , | ^| m |_>oo J- 

-K^ xl2 — -Ki xl2 -Ki ^2 — -*M ^2 

Since /i 2+ ^(8) is an intermidiate space, 

h 2+p (§) = (h 1+a (S),h 4+a (S)) {1+ ^ a)/3 , 

where (-|-) denotes the interpolation functor introduced by Da Prato and 
Grisvard [11], we obtain form [19, Proposition 2.4.1] that the elements on 
the diagonal of d<&(0) generate analytic semigroups, that is 

-d^i G H(h 4+f3 (S), h 1+a (S)),i = 1, 2, 

while the elements on the secondary diagonal belong to £(/i 2+ ^(§), /i 1+ ^ (§>)), 
that is they have lower order. We obtain then form [1, Theorem 1.6.1] that 
the matrix <9<£(0) is a generator, which completes the proof. □ 

We give now a short proof of our second main result, Theorem 1.2. 

Proof of Theorem 1.2. Let < j3 < a. Since <£ is analytic and <9<I>(0) gener- 
ates a strongly continuous and analytic semigroup, we find an open neighbor- 
hood O of in /i*+0(S) such that -d$( Pl ,p 2 ) G U({h i+ P(S)) 2 , (/i 1+ ^(§)) 2 ) 
for all (pi,p 2 ) G O 2 . Letting O := On /i 4+a (§)), we find that -d$(p 1 ,p 2 ) G 
^((/i 4+a (S)) 2 , (/i 1+Q (S)) 2 ) is, for all (p X ,p 2 ) G O 2 , the realisation of the op- 
erator -d$(pi,p 2 ) G %((/i 4+/3 (S)) 2 ,(/i 1+/3 (§)) 2 ). Whence, the assumptions 
of [19, Theorem 8.4.1] are all satisfied and the desired assertion follows at 
once. □ 
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